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INTRODUCTION 

Thus far, all analysis of system dynamics and the perfomance of process control systems 
have been carried out in either the time or Laplace domains. L^lace Transforms are used to 
expedite analysis of dynamical systems. Familiarity with Laplace Transform techniques 
means that we do not need to solve the underlying ordinary differential equations. The use of 
Laplace Transforms has been very usefiil and in many cases, analysis within the Laplace 
domain is more tban satis&^tory. However, there are other cases where this ^proach is 
inflexible and leads to intractable problems. In fact, frequency response techniques &ai be 
very useful in a whole spectrum of process control activities, e.g. 

• system identification 

• controller tuning 

• stability analysis 

• robusbie^ smalysis 

. design of noise filters 

SYSTEM roENTmCAHON 

After studying the fi^uency resp(»vses of various systems, it should be clear that the 

characteristics of the responses are determined by parameters of the systems. Not only do the 
sh^s and features of frequency response plots enable quick categorisation of system types 
and structures, they also enable more accurate determination of system parameters. Take the 
Bode diagram shown in Fig. 1 as an example. From the AR plot, we observe that there are 3 
asymptotes, one LFA and two HFA. The largest gradient of the high frequency asymptote 
(HFA2) is -40dB/decade and these combiae to indicate that the system is a second order 
system. The gain of the system is Kp and because there is no overshoot in the AR plot, it is 
most likely that the system is an overdamped one. Since it is an overdamped second order 
system, it should have two time constants, and these are given by the inverse of the 
fiiequencies at which the asymptotes intersect, i.e. Xi and T2. Now, if the system is purely 
second order, then the maximum phase shift should be -180°. However, the phase plot shows 
that the system's phase shift tends towards infinity, which is an indication that the system 
possesses a time delay. To determine the value of this time delay, all we need to do is to 
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Y(s) is the output; R{s) is the set-point; Kc is the proportional gain; Kp is process gain; the 
time delay, and x the time constant. Simplification gives 

Y(s) K^K^exp{-Qs) 



R(s) 1+15+ K^K^ exp(-e*) 

Tbe closed loop characteristic equation is therefore: 
1 + xs + K^K^ exp(-e5) = 

and the roots of this equation determines whether the closed loop is stable. The presence of 
the delay term therefore poses a difficulty, as it has to be approximated by a polynomial, to 
enable solution for the roots. Even discounting the accuracy of approximation, all that the 
vx^ tell us is whether the closed loop is stable or otherwise (absolute stability), not how near 
the system is to booming unstable (lektive stability). 

Frequency Domain 

By working in the frequency domain, the problem imposed by time delays in the transfer 
function ^)proach is avoided. Stability analysis in the fiiequency domain also tells us how 
close the system is to instability. Finally, the fi-equency domain allows the stability of the 
dc^ed loop system to be evaluated from the fi^quency response of the open loop system. Let 
us first exanune this last aspect. 

Let the Laplace transfer functions of the controller and the process be G^s) and Gp{s) 
respectively. The closed loq) duffact^l^ fcpiation is therefore: 

Any root that satisfy this equaticm also has to satisfy 

GXs)G^{s) = -\ 

In paiticulac, for strict stabihiy, the closed loop poles, must be such that 
G«(*)G,(5)U^<-1 

Transforming this alternative form of the characteristic equation into the fiiequency domain, 
we get 

G.C/«)G,0(0) = -1 
Expresang this in polar form, we obtain 

G,0®)G,0©)=K0Q)G,0to)|^G.a«)G^0«)) = |-l|^-l 
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Consider the polar plot above. Suppose the distance between the origin and the point where 
die fiequency response intersects the negative real axis is a. The value of the gain that will 
bring this to the -1 point is therefore l/a (since any additional gain is multiplicative). Thus the 
gain margin is given by the inverse of the distance between the origin and the point where the 
ftequency response intersects the negative real axis. 

To measure the phase margin from a polar plot, a circle of unit radius, and centred at the 
origin is drawn. A line is then drawn from the origin to the point where this circle intersects 
&e ftequency response. The angle between i3m line and the negative real axis will be the 
amount of extra phase shift that the system can tolei;^ before it beccnnes closed loop 
unstable, i.e. the phase margin, 

If all the open loop components of a system is grouped as Gol, then the detenninMion of gain 
and phase margins from polar plots can be expressed mathematically as: 

G„=r^ when ZGo^^-lSO" 

(|)„=180°+ZG^ when |g^| = 1 

Thus for the system to be stable, the gain margin must be greater than 1 and the phase margin 
must be positive. This also means that if the fi:equency les^nse crosses the (-1,0) point, then 
the closed loop system will be unstable. 
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Summary 

In summary, frequency response provides a powerful stability analysis tool. For a closed loop 
system to be stable, the frequency response of its open loop components must have a positive 
{diase maigin and the gain margin must be greater llian 1 (0 dB). It should also be clear that 
stable delay free first and second order systems under proportional control will always remain 
stable in the closed loop. This fact can best be visualised using the polar plot: the frequency 
response will never cross the negative real axis. 

Offline Z-N tuning 

Controllers can also be tuned using frequency response techniques. Perhaps the best known 
method is ftte Ziegler-Nidiols off-line tuning appioacb. bstead of attempting to drive the 
process to constant amplitude oscillation under proportional control, the offline approach uses 
a model of the open loop process to determines its gain maigin. Recall that the gain maigin is 
the amount of extra gain that can be introduced into the system before the closed loop 
becomes unstable. Thus, we know that when we impose this additional gain on the system, 
the closed loop will become marginally unstable, i.e. respond with constant amplitude 
oscillation. Thus, the gain maigin is equivalent to the "ultimate gain" determined 
e}q)6rimentaUy. The maiginally stable system will oscillate with a frequency equal to the gain 
cross-over fi^ency. Ihus, we can also determine the '\iitimate period". 

Ultimate gain: = 

2% 

Ultimate p^iod: P, = — ^M^^elHllg =gfflncross-ov<^frequeaK^ 

CO (J 

which is the inverse of the gain cross-over frequency. Having determined the "ultimate gain" 
and "ultimate pehod", these can be substituted into the Z-N tuning formulae below to obtain 
die impropriate P, PI or PID settiDgs for the process. 
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Ziegler-Nidiols Tuning Formulae 
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1 . Output Input 
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Figure 5. Result of relay feedback on output signal 

The ultimate gain is determined fix>ni, K^=4d/ita and Ifae ultimate period is given by the 
period of the output limit cycle, P„. Once these parameters are found, they can be used in the 
Z-N fomiulae to get ^propriate P, PI or PID settings. 

An advantage of online frequency domain based tuning techniques is that we do not need to 
know the structure of the process, i.e. whether it is first or second order, and whether it 
possesses a time delay. Unlike the online Z-N tuning method though, the relay feedback 
tuning strategy will always maintain the process aroimd its nominal set-point. The scheme 
presented here is also the basis for many commercial auto-tuning controllers. Hiese 
commercial devices use sUghtly more sophisticated approaches to mitigate the effects of 
process noise and disturbances, and variations of the Z-N tuning fermulae. A dismssion of 
how they function is however, beyond the scope of this course. 

ROBUSTNESS ANALYSIS 

In the design of a process control system, linear models of the plant are often used. The 
advantages are that: 

• linear models can readily be developed fixim plant operating data 

m there is a weatth of knowledge on tibe design of linear ccntrol systems. 
However, process plants are usually very complex systems, i.e. 

• they are affected by many inputs and outputs, some measurable some not 

• they suffirfiom varying time delays 

• they possess high order-dynamics 

• they often exhibit non-linear behaviour. 
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• the closed loop performance is acceptable (e.g. is tbe overshoot reasonable? is the 
rise time fast enough?) 

To investigate stability, we noimally develop the closed loop e}q)ression, i.e. 



Y{s) 




where R(s) is the desired output level or set-point 


Ris)' 


' l+G,is)G^is) 





• Closed loop stability is ensured if the roots of the characteristic polynomial, 
( 1 + Gr, (s)G^ (s) , have negative real parts 

• Given that the closed loop will be stable, closed loop performance can be easily 
assessed by ^plying the controller to a simulation of the linear model. 

If we are not h^py with either of the above results, then we have to re-design the P+I 
controller and repeat the analysis. 

Suiqx)se you are satisfied with the results of both investigations .... it is now time to apply the 
controller to the real process ! 

• What are the fectors that will compromise the integrity of the closed loop? 

The ccmtroUer has been designed based on a linear (assumed 
accurate) model of the process to be controlled. 

If the model is inaccurate or if the characteristics of the process 
change, then the designed control system will no longer be 
appropriate. 

This is the problem of process-model mismatch 

• In the presence of process-model mismatch, how confidmt am you that the closed 
loop will be stable and will perform as required? 

• Can we accoimt for these factors in iSbe design of Ihe controller prior to 
implementation? 

RobiutneM Aaa^s in the Laplace Domain 

Woiidng within the L^lace domain, Hh^ lUe two possible w^ of tackling the problem of 
process-model mismatch in control systems design. 

• Evaluate the designed controller against a variety of process models via simulation 
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(^^+A/:Jexp(-[e±Ae» 

1+(t±Ax)* 



= 



In this case, we would have to solve for fbe roots of Ifae closed loop equation for 
jmkm fxmhm&amQi ±AK^, ±A0 and ±A% l l 

Moreover, we would have to ensure that the magnitudes of AK^ , A0 and At are 
physically compatible. 

Robustness Analysis in the Frequency Domain 

Therefore, while Laplace Transforms have simplified the analysis of linear systems in many 
ways, its use becomes unwieldy when we want to investigate system properties under 'real' 
conditions. The fiiequency domain offers a much simpler ^pioach because the parameters of 
a process in the frequency domain are condensed into two parameters, the amplitude ratio 
and the phase -shift. This is regardless of the order of the process and whether the process 
possesses dead-time. Thus, imlike analysis in the Laplace domain, analysis of the effects of 
process model mismatch in linear systems reduces to studying tte effects on tiiese two 
parameters over a range of fiequoicy mi3l&isi. 

Id tiie fi^i^uency domain, the fiequency response of the s^'stem is simply a locus of fiequency 

dependent points plotted on the complex plane. For example, the frequency response of the 
nominal model is plotted as the continuous line in the figure below. Process model mismatch 
would lead to a dififeient locus of poinls. 
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Figure 6. Describing model un<»rtainties in the frequency domain 
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Figure 7. Bandwidth and cut-oflf frequency 

From the plot, we can observe that the amplitude ration decreases with increasing frequency. 
In other words, this first order lag has the ability to attenuate signals in the high fijequency 
lange. The degree of attenuation is usi^y tttmuie by the '1)andwidth", which is defined as 
the fi^quency range over which the amplitude ratio plot remains above -3dB. The upper 
fijequency value of the bandwidth is known as the "cut-off frequency" and defines the 
fiequency above which the signal is considered to be negligible. Note that the cut-off 
frequency is identical to the "comer frequency" for first order lags. Thus, we can deduce that 
a first order lag with larger time constants will have a smaller bandwidth, hence better 
filtering properties. This is shown below in the amplitude ratio plots of 3 first order lags, 
with time constants 10, 20 and 30. Here we can observe cleady that the bandwidth becomes 
smaller wi& increasing time conaMmts. 
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Ihe amount of phase shift that is incurred. Clearly, fiiequency domain analysis has a 
significant role to play in helping to solve this problem. 
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